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Purpose: To describe and validate a semi-automated targeted sampling (SATS) method for quantifying optic nerve 
axons in a feline glaucoma model. 

Methods: Optic nerve cross sections were obtained from 15 cats, nine with mild to severe glaucoma and six with normal 
eyes. Optic nerves were dissected, fixed in paraformaldehyde and glutaraldehyde, and processed for light microscopy by 
resin embedding, sectioning, and staining of axon myelin sheaths with 1% p-phenylenediamine before axon quantifica- 
tion. Commercially available image analysis software was used as a semi-automated axon counting tool (SCT) and was 
first validated by comparison with a manual axon count (MAC). This counting tool was then used in a SATS method 
performed by three masked raters and in a semi-automated full count (SAFC) method performed by a single observer. 
Correlation was assessed between the SCT and MAC using a linear model and analysis of covariance (ANCOVA). Cor- 
relation between the SATS and SAFC methods was calculated and the bias, systematic errors, and variance component 
assessed. The intraclass correlation coefficient (ICC) was determined to establish inter-rater agreement. In addition, the 
time required to perform the SATS and SAFC methods was evaluated. 

Results: Correlation between the axon counts obtained by the SCT and MAC was strong (r = 0.9985). There was evidence 
of an overcounting of axons by the SCT compared to the MAC with a percentage error rate of 13.0% (95% confidence 
interval [CI] 11.0%, 15.1%). Both the correlation of SATS count (average per rater) to SAFC (r = 0.9891) and inter-rater 
agreement (ICC = 0.986) were high. The SATS method presented an overall positive counting error (p<0.001) when 
compared to the SAFC, consistent with a fixed percentage overestimation of 11.2% (95% CI 8.3%, 14.2%) of the full 
count. The average time required to quantify axons by the SATS method was 10.9 min, only 27% of that required to 
conduct the SAFC. 

Conclusions: Our data demonstrate that the SATS method provides a practical, rapid, and reliable means of estimating 
axon counts in the optic nerves of cats with glaucoma. 



Glaucoma is a leading cause of vision loss in the adult 
human population and although less prevalent, a devastating 
cause of vision loss in children. Loss of vision results from 
the death of retinal ganglion cells (RGCs) and the loss 
of their axons, which constitute the optic nerve. It is now 
widely accepted that the initial insult responsible for the 
characteristic axonopathy occurs within the lamina cribrosa 
region of the optic nerve head. We have established a colony 
of cats with a spontaneously arising, recessively inherited, 
primary congenital glaucoma (PCG) associated with a muta- 
tion in the latent transforming growth factor beta binding 
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protein 2 (LTBP2; OMIM 602091) gene [1,2]. These animals 
represent an authentic homolog of human PCG at the Glau- 
coma 3, primary congenital, D (GLC3D; OMIM 613086) 
locus [3,4]. To facilitate the practical application of this 
model in translational research, an efficient, accurate, and 
reproducible method of quantification of RGC and/or optic 
nerve axon loss is required. Most conventional methods for 
quantifying RGCs are time intensive and costly, necessitating 
complex immune -labeling protocols or even requiring prior 
surgical intervention to back label RGCs with Fluoro-Gold 
(Fluorochrome, Denver, CO) or l.l'-dioctadecyl-3,3,3',3'- 
tetramethylindocarbocyanine perchlorate (Dil), which would 
be impractical and undesirable in our large animal model [5]. 
Quantification of RGCs in the feline retina is complicated 
by the substantial population of small (<10 |im) neurons, 
identified as displaced amacrine cells, that are estimated to 
account for approximately 80% of the cell population in the 
ganglion cell layer of the cat [6-9]. The histomorphometry 
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of RGCs is further complicated in felines by the uneven 
distribution of RGCs throughout the fundus, with the greatest 
density in the area centralis [10,11]. We therefore elected to 
focus our efforts on the evaluation of optic nerve axon loss in 
this model. The most accurate method of quantifying optic 
nerve axons is the full axon count technique performed either 
manually using transmission electron microscopy (TEM) or 
semi-automatically by light microscopy. These techniques are 
extremely time consuming and not 100% accurate, requiring 
extrapolations to circumvent their inherent underestimation 
of axon count numbers [12]. 

Faster methods, such as random sampling of the optic 
nerve cross-sectional area [13] or semiquantitative grading 
[12], present unreliable results due to focal variability in the 
degree of nerve damage and subjectivity of grading schemes. 
In our feline model, as in human glaucoma, a regional pattern 
of axon loss is often observed. Therefore, we elected to 
modify a rapid method, initially developed for use in rodent 
models, that combines automated counting and targeted 
sampling of distinct regions of relatively uniform optic nerve 
damage [14]. The purpose of our study was to validate a semi- 
automated targeted sampling (SATS) method for quantifying 
optic nerve axons in cats. 

METHODS 

Subjects: We used fixed optic nerve tissues collected 
postmortem from cats in a research colony that had been 
established from a pedigree of cats with spontaneously 
occurring recessively inherited PCG. Samples were selected 
from animals ranging in age from 6 months to 6 years and 
representing a range of different stages in the progression of 
disease. Weekly intraocular pressure (IOP) data, as measured 
by rebound tonometry [15], were available for all cats in the 
study. All glaucomatous animals had persistently elevated 
IOP, a mean (standard deviation [SD]) IOP of 41.4 (14.6) 
mmHg, while mean (SD) IOP in normal cats was 18.3 (2.4) 
mmHg. The cup-to-disc ratio cannot be reliably assessed 
on ophthalmoscopy in cats, but all cats in the former group 
exhibited abnormalities consistent with glaucoma on elec- 
trophysiological testing (visual evoked potentials and pattern 
electroretinography) or optical coherence tomography. Single 
optic nerves were selected from normal cats (n = 6) and from 
cats that demonstrated mild to severe PCG due to a consistent 
mutation in LTBP2 (n = 9). All procedures were conducted 
with the approval of the University of Wisconsin-Madison 
institutional animal care and use committee. 

Tissue processing: Animal were euthanized by an intravenous 
overdose of pentobarbital sodium (120-180 mg/kg). Imme- 
diately following euthanasia, animals were transcardially 
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perfused with ice-cold 0.1 M PBS (100 mM phosphate, 154 
mM sodium chloride, pH 7.4), followed by 4% paraformalde- 
hyde solution in 0.1 M PBS (Electron Microscopy Sciences, 
Washington, PA). Eyes were enucleated and fixed overnight 
in 4% paraformaldehyde in 0.1 M PBS (Electron Microscopy 
Sciences) at 4 °C, then transferred to 0.4% paraformaldehyde 
in 0.1 M PBS for storage at 4 °C. An approximately 2-mm 
long segment of optic nerve was dissected 2 mm posterior 
to the globe. Shallow radial orientation cuts were made, one 
superiorly and two temporally, before postfixation in 2.5% 
glutaraldehyde/0.1 M PBS for 48 h at 4 °C. Nerve samples 
were then osmicated in 1% osmium tetroxide in 0.1 M 
PBS, rinsed, and dehydrated through an ascending series 
of alcohol treatments before routine epoxy resin embed- 
ding and sectioning. Semithin (l-|am) sections were stained 
with 1% p-phenylenediamine (PPD) for evaluation by light 
microscopy. 

Optic nerve imaging and axon counting: Full optic nerve 
cross sections were analyzed from each optic nerve. Images 
from the PPD-stained optic nerve sections were obtained 
using a bright light microscope Olympus BX43 (Olympus 
Inc., Center Valley, PA) with attached Olympus DP72 digital 
camera (Olympus Inc.) and captured using commercially 
available image analysis software (cellSens Dimension®, 
Olympus Inc.). As described below, the software's semi-auto- 
mated axon counting tool (SCT) was compared to a manual 
axon count (MAC) for validation and then used in SATS and 
semi-automated full count (SAFC) methods. Comparisons 
were then made between the SATS and SAFC methods for 
axon counting. All counts were performed by raters masked 
to animal identity and disease status and to the other raters' 
results. 

Validation of the semi-automated counting tool: A semi- 
automated axon detection process using the built-in auto- 
matic object recognition tool of the image analysis software 
was compared to the MAC method by a single rater. Digital 
images from nine different fields were obtained under a 63X 
microscope objective lens. The images represented three 
subjectively distinct areas of the optic nerves from each of 
three different affected animals including those judged to 
be normal or densely populated areas (n = 3), moderately 
damaged areas (n = 3), and severely damaged areas (n = 
3). Semi-automated axon detection was performed using 
cellSens Dimension® image analysis software. The images 
were analyzed under the "count and measure" menu. Using 
the adaptive threshold tool, the black PPD-stained myelin 
sheath of the axons was selected by visually adjusting the 
threshold values until all the axons in the field were selected, 
taking care not to obliterate the smaller axons (Figure 1). 
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We established that a threshold with a minimum intensity 
value of -255 and maximum value of -11 was appropriate for 
our samples, and this threshold was used for all subsequent 
semi-automated axon counting (including SATS and SAFC 
methods described below). The "auto split objects" tool was 
applied to the whole image to separate contiguous axons, and 
the "count and measure" tool was used to count the axons. 
Nonaxonal darkly staining elements of the optic nerve tissue 
that were being counted by the software were manually 
selected and deleted. Axons at the margins of the image that 
were not entirely visible were automatically excluded from 
counting by the software. The final axon count for each 
image was recorded. 

Manual axon counts: Manual axon counts were obtained for 
the same nine digital images by the same rater. Each and 
every axon in the image was visually identified and selected 
using the software's manual count tool. The results obtained 
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by the SCT were then compared to the MAC to determine the 
agreement between the techniques. 

Semi-automated targeted sampling method: The SATS 
method was adapted from the method previously described 
by Marina et al. [14]. Briefly, up to three distinct relatively 
homogeneous regions were visually identified in low-power 
(10X objective) photomicrographs of the entire optic nerve 
cross sections, and the area of each region (in um 2 ) was 
recorded using the image analysis software (Figure 2). 

Axon counts from up to five randomly selected high- 
power fields (63X objective) from each region were obtained 
using the SCT described earlier. The average axon density 
(axons/um 2 ) for each region was calculated by dividing the 
average axon count of the selected fields by the known high- 
power field area (242,611 um 2 ). Then, by multiplying the 
average axon density of that region by the total region area 




t 



4 ;* 







fx 



,3* 



v.OJ 



Figure 1. Semi-automated counting 
tool used in axon detection and 
quantification. Histologic sections 
of the optic nerves from affected 
cats. The images represent three 
distinct areas of the optic nerves 
judged as normal or densely popu- 
lated (A), moderately damaged 
(C), and severely damaged (E). 
Images B, D and F represent the 
corresponding images after semi- 
automated axon detection was 
performed using the image analysis 
software. The contour of the black- 
stained axon myelin sheath is high- 
lighted in red and automatically 
counted by the software. The scale 
bar = 20 urn. 
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Figure 2. Selections of regions 
with homogeneous damage in an 
optic nerve cross section from a 
glaucomatous cat. Images A, B, 
and C depict region selection by 
three different raters in the same 
optic nerve section. The blue region 
labeled (1) represents a normal 
region, the green region labeled (2) 



a moderately damaged region, and the red region labeled (3) a severely damaged region. Note the variation of the selected areas by the 
different raters and the final area of each region in um 2 . The scale bar = 100 um. 



((j.m 2 ), we obtained the estimated axon count for each region. 
Estimated counts for each distinct region were subsequently 
summed, and an estimate of the total count over the entire 
cross section was thus obtained. 

To estimate inter-rater variability, three different masked 
raters applied the SATS method to the same set of 15 optic 
nerve cross sections (nine affected cats and six controls). 
Raters selected regions and representative fields within 
regions independently. In addition, one rater recorded the 
time required for image acquisition and counting using both 
the SATS and the SAFC methods (see below). 

Semi-automated full count method: SAFCs were obtained 
by a single masked rater for the same set of 15 optic nerve 
cross sections in which the SATS method was applied. A 
consecutive series of images (up to 40 images) representing 
the entire optic nerve cross section was captured under a 63X 
microscope objective. To ensure the whole tissue section was 
captured without duplications or omissions, anatomic and 
other tissue landmarks, such as blood vessels or indentations, 
were used as reference points while manually navigating 
across each slide. All axons in each image were then counted 
using the SCT, as described previously, and the sum of the 
counts from all photomicrographs provided a total count of 
axons for the entire nerve cross section. The time taken to 
acquire images and axon counts by this method was recorded. 

Statistical analyses: To validate the SCT, we calculated the 
correlation between the axon counts obtained by the SCT and 
by the MAC techniques in the nine different images and we 
fit a linear model of the counting error (SCT less MAC) to the 
MAC. An analysis of covariance (ANCOVA) was performed 
to determine if this linear relationship varied by slide, and we 
performed tests for non-zero intercept and non-zero slope. 

To compare the axon counts obtained by the SATS 
method performed by three raters for the 15 optic nerves to 
the SAFC for the same nerves, we fit a linear mixed-effects 
model of the counting error (SATS count less SAFC count) 



that included a fixed-effect linear relationship to the SAFC 
count and random effects for slide and rater. We tested the 
SATS method for an overall bias using a Wilcoxon signed 
rank test of the differences between the average SATS count 
(averaged over raters) and the SAFC count for each slide, and 
we fit a two-stage linear model of the average counting error 
(average SATS count less SAFC) versus SAFC to test for non- 
zero slope and non-zero intercept and to estimate constant 
absolute error and constant percentage error models for the 
counting error. We also used the linear mixed-effects model 
to calculate variance component estimates. We determined an 
ICC for inter-rater agreement, and we estimated the correla- 
tion of the SATS method to the SAFC method for a typical 
rater by averaging the correlations between SAFC and SATS 
counts for each rater. 

We also performed a comparison of the time required 
for one reviewer to count with the SATS and SAFC methods 
using a Wilcoxon signed rank test. Data analyses were 
performed using R version 3.0.1 [16] and the lme4 package 
version 1.0-4 [17]. 

RESULTS 

Comparison of counts of semi-automated axon counting tool 
to manual axon counts: Figure 3 shows a scatter plot of axon 
counts obtained by SCT versus MAC for the nine fields with 
the least-squares linear fit (solid) and a line representing 
equality of the counts (dotted) for comparison. The correla- 
tion between SCT and MAC was high (r = 0.9985). 

An ANCOVA analysis showed no significant differences 
by slide in the fitted intercept or slope. Also, there was no 
evidence of a non-zero intercept, consistent with a constant 
percentage error rate over the large range of per- field counts 
sampled (614 to 2,483 axons). A model with zero intercept 
gave an estimated constant percentage error rate of 13.0% 
(95% CI 11.0%, 15.1%), representing overcounting of axons 
by the SCT compared to the MAC technique. Individual raters 



379 



Molecular Vision 2014; 20:376-385 <http://www.molvis.org/molvis/v20/376> 



© 2014 Molecular Vision 



CO 
O 
O 
O 



O 

CO 

o 
o 



3 

o 
O 

c 
o 
x 
< 

T3 
CD 

-t— > 

CO 

E 
o 
-*— ' 

Zi 
CC 
J_ 

E 

CD 

CO 



CO 

c O 
CD o 
CD O 

-■ 

CD ,., 
CD ^ 

£ -a 

CD CO 

Is 

CO 



o 

CO 



LO 
C\i 



o 
c\i 



o 



q 

O 




Manual Axon Count (MAC) in 1000s 



co 




Average of SCT and MAC (1000s) 



Figure 3. Degree of agreement 
between semi-automated axon 
count tool (SCT) and manual axon 
count (MAC) methods. The top 
panel shows a scatter plot of SCT 
versus MAC counts for n=9 fields, 
one each of low, medium, and 
high density selected from each of 
three animals (designated by three 
different plotting symbols). The 
dotted line represents equality of 
counting methods, and the solid 
line represents the best linear fit. 
The bottom panel shows a Bland- 
Altman plot of the difference versus 
the average of the two counting 
methods. 



selected areas and representative fields independently and 
subjectively, and although choice of areas and fields varied 
considerably from rater to rater (Figure 2), both inter-rater 
agreement (ICC = 0.986) and the average per-rater correlation 
of SATS count to SAFC (r = 0.9891) were high. 



The mean of the SAFC for the 15 optic nerves sampled 
was 52,827 axons with an SD of 24,483 axons. Figure 4 shows 
a scatter plot of SATS counts by three raters for each of the 
15 optic nerves versus the SAFC for those nerves performed 
by one of the raters. A solid line giving the best linear fit is 
shown as is a dotted line representing equality of the counts 
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Figure 4. Degree of agreement 
between semi-automated targeted 
sampling (SATS) and semi-auto- 
mated full count (SAFC) methods. 
The top panel shows a scatter plot 
of SATS counts performed by 
three raters versus the SAFC count 
for each of n=15 optic nerve cross 
sections. Raters are identified by 
plot symbol color, the dotted line 
represents equality of the counting 
methods, and the solid line repre- 
sents the best linear fit. The bottom 
panel shows a Bland-Altman plot of 
the differences versus the average 
of the two counting methods. 



for comparison (line of equivalence). The fitted linear rela- 
tionship is 

SATScount = 1 5 8 + 1 . 1 098 x SAFC 
or equivalently: 



SA TC - counting - error = 158 + 0.1 098 x SAFC 



representing strong evidence of an overall positive 
counting error (p<0.001) with overcounting of the SATS 
method compared to the SAFC method. There was no 
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evidence of a non-zero intercept, consistent with a fixed 
percentage error rate of 11.2% of the total count (95% CI 
8.3%, 14.2%). 

In the fitted linear mixed-effects model, the random 
rater effect had zero variance, indicating consistency of 
the observed error pattern across raters (i.e., no systematic 
tendency of a rater to over- or undercount axons compared 
to other raters). The estimated residual error SD was 3,287 
axons per slide (95% CI 2,603; 4,333), which represents the 
rater-to-rater variation observed in SATS counts performed 
by different raters of a single optic nerve slide. The estimated 
SD of the slide random effect was 2,572 axons per slide (95% 
CI 862; 4,054), which represents additional slide-to-slide 
variation attributable to features of a particular optic nerve 
cross section that make it more subject to over- or under- 
counting by the SATS method compared to other sections 
with similar axon density. 

Comparison of time spent performing semi- automated 
targeted sampling and semi-automated full count methods: 
For the rater who applied both the SATS and SAFC methods 
to the 15 optic nerves, the SATS method took significantly 
less time, on average only 27% as much time, as the SAFC 
method: SATS mean (SD) 10.9 (4.9) min versus SAFC mean 
45.3 (10.7) min; (p<0.001, n=15, paired Wilcoxon signed rank 
test). This included time saved in both image acquisition 
(mean 6.3 versus 21.1 min) and actual counting time (mean 
5.7 versus 24.3 min). 

DISCUSSION 

To validate the SCT, we analyzed nine fields of optic nerve 
cross sections with variable density and compared the 
obtained axon counts to MACs of the same fields. The corre- 
lation between the two approaches was high (r = 0.9985), 
but there was a consistent bias — an overestimation of axons 
by 13% when using the SCT compared to the MAC. The 
most likely explanation is a misidentification of axon-like 
structures by the software's automatic detection system. 
On any given optic nerve cross section, and particularly in 
moderate to severely affected areas, myelin or cellular debris 
from degenerating axons are stained black by PPD and can 
be recognized and counted by the software. Furthermore, 
splitting of myelin fibers can cause a discontinuity of the 
circumferential myelin sheath surrounding an axon, causing 
the software to recognize the axon as two separate objects. To 
circumvent these problems we established an additional step 
in the protocol that includes manually selecting and deleting 
the detected nonaxonal structures and double-counted axons. 
Despite this extra effort, some of the artifacts erroneously 
counted by the software could not be easily recognized and 
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excluded by the rater. The addition of this step caused a slight 
increase in the overall time spent analyzing a slide, but even 
with its flaws it proved to be important in creating a homo- 
geneous set of data, a fact that can be appreciated by the high 
correlation between the SCT and MAC. 

Our ANCOVA analysis showed that this 13% overesti- 
mation did not vary significantly across the three slides (that 
is, there was no tendency for any slides to be more or less 
prone to overestimation than the others). This suggests that it 
would be possible to develop an accurate calibration curve to 
correct for this bias in the SCT. Such a correction would need 
to be developed separately for a particular slide preparation 
and counting protocol using a larger sample of slides than we 
used in our analysis. However, our ANCOVA results provided 
confidence that even uncorrected counts obtained by the SCT 
are comparable from slide to slide because of the consistency 
of counting bias. 

Full optic nerve axon counts in normal cats have been 
reported in the literature with values ranging from 86,000 
axons by light microscopy [18] to 193,000 axons by TEM [19]. 
Control cats from our study presented mean axon counts of 
73,751 axons with the SAFC method and 78,380 for the SATS 
method (average of three raters). It has been reported that 
light microscopy techniques tend to underestimate total axon 
counts in the order of 20%-30% [14]. This could be explained 
by the capacity of TEM to detect small axons that were most 
likely not detected by the SCT. Approximately 48%-52% 
of cat axons present diameters between 0.3 to 1.5 um [20], 
and those with diameters in the lower end on this spectrum 
might not be detected by our method. In addition, the mesh 
grid employed in TEM procedures occupies 55% of the image 
area, leaving only 45% of the tissue available for an axon 
count; thus the total axon count must be extrapolated, with 
the potential for overestimation of axon counts. The cat optic 
nerve presents marked fascicular organization of axon fibers 
[21] as well as size grouping and regional variations in axonal 
density [22]. These factors could impact the final axon counts 
using TEM since the extrapolations assume that the remaining 
55% of the tissue is of similar composition, introducing the 
potential for overestimation of axon numbers. However, it 
was not the intent of our study to determine reference ranges 
for axon counts in cats but to develop a reproducible and reli- 
able method to quantify axons in our model. These counts, in 
turn, serve as a basis for validation in individual cats of the 
noninvasive structural and functional measures of glauco- 
matous damage (acquired by optical coherence tomography 
and electrophysiology, respectively) that are the focus of our 
ongoing longitudinal studies. Analysis of the optic nerves of 
a much larger number of subjects is now underway and might 
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allow us to better characterize a reference range for feline 
axons counts obtained by light microscopy. 

Both the correlation between the SATS and SAFC 
methods (r = 0.9891) and the inter-rater agreement (ICC = 
0.986) were high, confirming the accuracy, validity, and 
reproducibility of the SATS method in our animal model. 

The SATS method presented an overall positive counting 
error when compared to the SAFC method, consistent with 
a fixed percentage overestimation (11.2% of the full count). 
One explanation for the overestimation of axon counts using 
the SATS method is that the mean axon count of the five 
selected fields is multiplied by the region's area (urn 2 ). This 
calculation assumes that each field is composed only of 
axons. In reality, optic nerve sections of normal cats included 
in this study demonstrated that approximately 15% of the area 
of a high power (63X) field is composed of nonaxonal tissue 
(data not shown), and thus overestimation of axon numbers 
is to be expected. This problem is hard to address since to 
correct for this source of overestimation the nonaxonal area 
needs to be extracted from the images of each selected field 
and region. This process is time consuming and would defeat 
our initial purpose of creating a more rapid efficient method. 
Another alternative would be to apply a correction factor to 
the final counts. However, this correction factor would be 
hard to establish given the heterogeneity in different optic 
nerve regions and the fact that in the most affected areas it 
is hard to distinguish areas that represent optic nerve native 
connective tissue from those that represent areas with genuine 
axon loss. 

When using time-efficient approaches, like semiquanti- 
tative grading schemes [12,23,24], random sampling [25,26], 
targeted [14,27], or automated full counts [28], to assess 
optic nerve damage, the subjectivity of important steps in 
the protocol can introduce significant variability in the data. 
Many of the published studies using these counting methods 
and grading schemes have not addressed variance between 
raters and optic nerve slides or attempted to evaluate sources 
of bias in counting. In our study, raters selected homogeneous 
optic nerve regions and fields independently and subjectively, 
and there was marked variability in the selections made by 
our raters (Figure 2). However, despite this subjective step, the 
inter-rater agreement (ICC = 0.986) and correlation of axon 
counts obtained by SATS (average per rater) to those obtained 
by SAFC (r = 0.9891) were high. The data analysis also found 
no significant random rater effect, meaning no tendency for a 
given rater to systematically over- or undercount axons. This 
suggests that, assuming adequate training in its application, 
the SATS method might be usefully employed in a large 
study extending over several years where it is not feasible 



© 2014 Molecular Vision 

for a single rater (or the exact same set of multiple raters) to 
perform every SATS count. 

In experiments that generate a large number of optic 
nerves for analysis, time becomes an important consideration. 
Manual counting of axons on TEM and light microscopy 
images can be time consuming to the point of becoming 
impractical. In the current study, the average time spent 
acquiring images, selecting regions and fields, and counting 
axons for the SATS method was less than 30% of that required 
to conduct SAFC. These time savings were facilitated by the 
commercially available image analysis system used that, 
after appropriate thresholds were set, delivered axon counts 
in seconds. Alternative image analysis software packages, 
whether commercially or freely available, present similar 
capabilities and could be used with appropriate adaptations 
and validation. 

Our method has several limitations. First, as for any other 
image analysis protocol, the results depend largely on the 
quality of the images analyzed. Image quality is influenced by 
proper processing, sectioning, mounting, and staining of the 
optic nerve tissue. Second, while we have provided character- 
izations of and models for the observed counting error when 
comparing SCT with MAC and SATS with SAFC methods, 
these errors were observed under a specific and strict slide 
preparation and analysis protocol that included, among other 
things, selection of specific values for the software's detection 
"threshold" for the SCT. Counting errors are almost certainly 
dependent on details of the counting protocol, including the 
selected threshold value, and so our error models must be 
considered specific to the protocol and detection settings 
used. In particular, our models of counting error are unlikely 
to be useful to "recalibrate" SATS counts to match manual 
counts in a more general setting. Third, because our results 
were based on a single set of slides prepared and counted 
over a relatively short period of time according to a specific 
strict protocol (including specific threshold value), we were 
not able to fully characterize errors that might be observed 
across separate slide preparations at different time points, as 
might occur in a large multiyear study. With respect to the 
SCT and MAC comparison, while the correlation between the 
two methods was high and we found no evidence of slide to 
slide variation in the error rate, the sample number was small, 
encompassing only three animals and nine fields. 

In summary, our data indicated that the SATS method 
provides a practical, rapid, and reliable means of estimating 
axon counts in the optic nerves of normal cats and cats with 
glaucoma. Despite the relative subjectivity of certain aspects 
of the method, it provided reproducible axon counts that were 
comparable to those obtained by the SAFC method while 
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reducing the time taken to acquire axon counts by more than 
two-thirds. 
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